Asymptotic Capture-Number and Island-Size Distributions 
for One-Dimensional Irreversible Submonolayer Growth 
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Using a set of evolution equations [J.G. Amar et al, Phys. Rev. Lett. 86, 3092 (2001)] for the 
average gap-size between islands, we calculate analytically the asymptotic scaled capture-number 
distribution (CND) for one-dimensional irreversible submonolayer growth of point islands. The 
predicted asymptotic CND is in reasonably good agreement with kinetic Monte-Carlo (KMC) results 
and leads to a non-divergent asymptotic scaled island-size distribution (ISD). We then show that a 
slight modification of our analytical form leads to an analytic expression for the asymptotic CND 
and a resulting asymptotic ISD which are in excellent agreement with KMC simulations. We also 
show that in the asymptotic limit the self-averaging property of the capture zones holds exactly 
while the asymptotic scaled gap distribution is equal to the scaled CND. 

PACS numbers: 68.55.Ac,68.43. Jk,89.75.Da 
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Recently, considerable theoretical effort has been car- 
ried out towards a better understanding of the scaling 
properties of the island-size distribution in submono- 
layer epitaxial growth iiS42iiiS4S4L& For example, in the 
pre-coalescence regime the island size distribution Ns{6) 
(where Ng is the number of islands of size s at coverage 
9) satisfies the scaling form,^ 
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where S is the average island size, and the scaling func- 
tion f{u) depends on the critical island size and on the 
island morphology.^ 

One of the standard tools used in these studies is the 
rate-equation (RE) approach which involves a set of de- 
terministic coupled reaction-diffusion equations describ- 
ing the coverage-dependence of Ns{9) through a set of 
rate-coefficients usually called capture numbersAiS For 
the irreversible growth of point islands, rate-equations 
valid in the pre-coalescence regime may be written in the 
form 
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where the capture numbers a, (ui) correspond to the 
average capture rate of diffusing monomers by islands 
of size s (monomers) and R = D/F is the ratio of the 
monomer diffusion rate to the deposition rate. Accord- 
ingly, the central problem in using the RE approach is the 
determination of the average capture numbers a, (9) and 
the corresponding capture number distribution (CND). 

Recently, we have developed a self-consistent rate- 
equation approach to irreversible submonolayer growth 



in which correlations between the size of an island and 
the corresponding average capture zone are explicitly 
taken into account in order to accurately predict the 
scaled island size and capture number distributions. 
Our method involves numerical integration of the island- 
density RE's (2) along with the analytical solution of 
a set of approximate evolution equations for the average 
Voronoi-area surrounding an island of size s, and is based 
on the following two assumptions: (1) the average cap- 
ture area per freshly nucleated dimer is proportional (be- 
fore rescaling) to the average area per island and (2) the 
combined effects of the preferred nucleation of dimers in 
large capture zones and the preferred "break-up" of large 
capture zones due to nucleation may be approximated by 
a uniform rescaling of the average capture zone of each 
island&i Using this approach, we have obtained numeri- 
cal results for the scaled capture number and island-size 
distributions which agree well with kinetic Monte Carlo 
(KMC) simulations in both one and two dimensions over 
a wide range of experimentally relevant values of D/F 
{D/F = 10^ - 10^). 

Recently, it has been argued™ that because our 
method docs not explicitly take into account spatial fluc- 
tuations in the nucleation, then it must lead to a diverg- 
ing ISD in the asymptotic limit corresponding to infinite 
D/F. However, as shown by Bartelt and Evans, ^ in the 
asymptotic limit the scaled ISD is related to the scaled 
CND as, 



f{u) - /(O) exp 
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where C{s/S) — Gg/oav is the scaled CND, z is the dy- 
namical exponent describing the dependence of the aver- 
age island size on coverage {S ~ 9^), and /(O) is deter- 
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mined by the normalization condition, 

/•oo 

/ du f{u) = 1. (4) 

As pointed out in Ref. 2, Eq. Elimplies that if C(u) > zu 
then no divergence will occur. However, if C{u) crosses 
zu at some value Uc then the ISD will be cut off at Uc (i-e. 
f{u) = for u > Uc) if C"(zic) > 2z — 1 while a divergence 
in the ISD will occur if C'{uc) < 2z — 1. An example of a 
divergent asymptotic ISD is the usual mean-field theory 
with C{u) = 1. Thus the question of whether or not our 
method leads to a divergence in the asymptotic limit is 
entirely determined by the asymptotic scaled CND. 

In this Rapid Communication we rigorously address 
the question of the asymptotic behavior obtained us- 
ing our method by analytically deriving the asymptotic 
scaled CND along with the corresponding asymptotic 
scaled ISD for the case of irreversible growth of point 
islands in one-dimension. We find that, contrary to the 
claims in Ref. 11, our method leads to a non-divergent, 
asymptotic ISD as well as to an asymptotic scaled CND 
which is close to that obtained in simulations. We then 
show that by slightly modifying our original analytical 
form for C (u) , an improved analytical expression for the 
asymptotic scaled CND may be obtained which is in ex- 
cellent agreement with KMC simulations. The resulting 
analytic expression leads to a scaled ISD which is also 
in excellent agreement with KMC simulations. We also 
demonstrate that the asymptotic scaled gap distribution 
is identical to the scaled CND and as a result the self- 
averaging property of the capture zones holds exactly. 

For clarity, we briefly review our method and its appli- 
cation to the case of irreversible growth of point islands 
in one-dimensioni In this case, we have shown that the 
local capture- number (7(j/) for an island with gap-size y 
(corresponding to the distance to the nearest island) is 

a(y;0) = ^ (4) ' 

where the monomer capture length ^ and the nucleation 
length ^1 are defined as 

2ai7Vi = 1/C?, 1/e? + ^ (7, TV, = 1/e. (6) 

From the evolution equations for the distributions of gap 
lengths, the average "gap-length" j/s (before rescaling) 
corresponding to an island of size s may be obtained as 
the solution of the equation, 

s-2 = / i?iVi(0)a(y,;</))d(/). (7) 

The coverage 9y is defined by = 6 Y{9y) where b is the 
proportionality factor (before rescaling) between the av- 
erage gap-length of a freshly nucleated dimer and the av- 
erage gap-length of all islands, Y{d) = 1/N{9) is the av- 
erage gap-length at coverage 9, and N{9) = X]s>2 ^s{9) 



is the average island density. Physically, the integral in 
Eq. [71 may be interpreted as corresponding to the aver- 
age number of particles added to the dimer since it was 
formed at coverage 9y , neglecting any change in the cap- 
ture zone due to nucleation and break-up. To include 
the effects of break-up, the gap lengths are rescaled by a 
rescaling factor a — l/X]s>2-^s Vs so that the average 
gap length for an island of^ size s is given by = a ys 
while the capture number is given by (t, = (j{ys)- 

We now consider the asymptotic limit corresponding to 
infinite D/F. In this limit, and assuming that 9 » 9^ 
(where 9x ~ R~^^^ corresponds to the coverage at which 
the island-density equals the monomer density), Eq. [S] 
may be rewritten as, 

aiy;9)^y/e. (8) 

This implies that in the asymptotic limit the scaled gap- 
distribution B{s/S) — ys/Y is identical to the corre- 
sponding scaled capture number distribution C{s/S). 
This result also demonstrates the "self-averaging prop- 
erty" of the capture zones, i.e. in the asymptotic limit 
the average capture number as is exactly equal to the 
local capture number evaluated at the average capture 
zone for an island of size s. 

In the asymptotic limit that 9 >> 9x one also has 
dNi/d9 = l-RNx/^"^ ~ 0, which implies that RNi/(^ = 
1. Using these results in Eq. [7|and assuming that also 
9y » 9xr^ we obtain 

s^2^ys{9~ 9y). (9) 

Since in the asymptotic limit 5^1, this may be rewrit- 
ten as 

u = s/S^{ys9/S){l~9yl9). (10) 

In the asymptotic limit, one has 9/ S ^ N ~ ^/Y. Us- 
ing = b/N{9y) and N{9y) ~ 91^'^, N{9) - 9^/^ (since 
9y, 9 ^ 9x and z = 3/4: for irreversible growth of point- 
islands in one-dimension), one has {9y/9) = {bY/ijs)'^. 
Thus, the asymptotic gap-distribution before rescaling 
will satisfy, 

u = B{u) -b^/B{uf, (11) 

where B{u) — ys/Y. Rescaling, we obtain for the asymp- 
totic scaled gap distribution B{u) — a B[u) the equation 

u = B{u)/a - b"- [B{u)/a]-^, (12) 

where the rescaling factor a and the proportionality fac- 
tor b must satisfy the requirement that 

poo 

/ B{u)f{u) du = 1 (13) 

It is easy to see that for any positive constants a and 6, 
Eq. 12 has a unique positive solution for every positive 



u while for a > z one has B{u) > zu. Similarly, one may 
show that for a < 2/3 the ISD predicted by Eq. 3 will 
diverge while for 2/3 < a < z there will be a cut-off at 
Uc, i.e. /(m > Uc) = 0. 

We now consider the case = 1 as was assumed in our 
numerical calculations for finite D/F.^-"^ Eq. |2|along with 
the requirement that f{u) — > rapidly for large m, indi- 
cates that the asymptotic B{u) should satisfy B{u) ~ zu 
for large u& Applied to Eq. ^] this implies that a — z 
and thus C(0) — B{0) = 3/4, in very good agreement 
with our KMC simulation results. This leads to the fol- 
lowing expression, 

u = I3{u)-1/B{uf (14) 

where B{u) = B{u)/z and B{Q) = 1. Unfortunately, 
integrating Eq. 13 numerically using Eq. 1141 to deter- 
mine f{u), we find that the resulting "cross- integral" 
B{u) f{u) du ~ 1.15 is somewhat larger than 1. This 
implies that a smaller value of the rescaling factor a must 
be used. By carrying out the appropriate numerical inte- 
grals we find that a ~ 0.70 satisfies the "cross-integral" 
normalization condition H13|) . The corresponding scaled 
CND and ISD with 6=1 are shown in Fig. 1. As can 
be seen, the calculated asymptotic ISD is shifted to the 
right from the simulation result and is cut off to zero at 
Uc — 1.85 where C(u) crosses zu. Thus we find that, 
contrary to the claims in Ref. jO, our method leads to 
a non-divergent ISD in the asymptotic limit of infinite 
D/F. 

It is also interesting to consider the general case b 1, 
i.e. the capture-area of a freshly-nucleated dimer is pro- 
portional to the average area per island with some un- 
known proportionality constant b. In this case, we as- 
sume that the value of the rescaling constant a is fixed 
to the value a = z by the requirement that the asymp- 
totic behavior of B{u) is given by B{u) ~ zu for large 
u. We then search for a value of 6 (6 ~ 0.87) such that 
the sum-rule (|13|) is satisfied. As can be seen in Fig. 
1, the resulting scaled gap-distribution B{u) is now very 
close to the simulation results for all u except for small 
u {u < 0.7) where it is now significantly lower than for 
the case b—1. Accordingly, the peak of the correspond- 
ing island-size distribution (Fig. 1(b)) is somewhat lower 
than the peak of the simulated distribution, and the dis- 
tribution itself is somewhat wider. 

We now present an analytical form for the capture- 
number distribution C{u) — B{u) which provides excel- 
lent agreement with simulations and thus strongly sup- 
ports our conjecture that -B(O) = z and that B{u) zu 
for large u. Since Eq. E|with b =1 and a ^ z leads to 
an expression for B{u) H14() which satisfies B{0) = 1 and 
B{u) zu for large u, but a cross-integral H13|l which 
is just slightly larger than 1, we propose that the cor- 
rect analytic expression for B(u) has the same form as 
(|14ll but with additional terms corresponding to higher 
order powers of [B{u)]~^. As an example, we chose the 
following form, 

u^ B{u)-{b2/B{uf) exp[&i/B(u)"], (15) 
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FIG. 1: Calculated asymptotic scaled (a) gap-distribution 
B{u) and (b) island-size distribution f{u), along with KMC 
results (symbols). 



where 62 = e^''^ in order to satisfy B{0) = 1, and n is a 
free parameter for which we will choose n — 5. We then 
use Eq. |31to calculate f{u) and search for a value of 61 
(&i ~ 0.993) such that the sum-rule ((T^ is satisfied. Our 
results for the corresponding B{u) and f{u) obtained in 
this manner are shown by the solid curves in Fig. 1. As 
can be seen, the agreement between the calculated scaled 
gap-distribution B{u) and simulations is now excellent. 
Accordingly, the agreement between the corresponding 
calculated ISD and simulations is also very good. This 
result also confirms that C{u) = B{u) in the asymptotic 
limit. 

We have also tried several other values of n (n = 
3, 4, 6) and obtained similarly good agreement with 
simulations.^"* Finite expansions in powers of [B{u)]~^ 
satisfying B{Q) = 1 also lead to reasonably good agree- 
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ment. In particular, the expression 

u = B{u) ~ lB{u)-^ - + B{u)-'% (16) 

(with 7 ~ 0.39 to satisfy H13|l ) leads to results which 
are essentially identical to those obtained using Eq. [Tf)l 
Thus, it would appear that almost any form similar to 
(|14|l . but with higher order corrections, which satisfies 
the condition B{Q) = 1 as well as the sum-rule (|13|l . 
leads to reasonably accurate results for the asymptotic 
scaled ISD and CND. 

In conclusion, we have derived analytical expressions 
for the asymptotic scaled capture-number, gap-, and 
island-size distributions for the case of one-dimensional 
irreversible growth of point islands. In particular, we 
have shown that our method leads to an asymptotic 
scaled CND which is close to that obtained in simula- 
tions and as a result to a non- divergent asymptotic ISD. 
Our analytical results also indicate that in the asymp- 
totic limit the scaled gap-distribution B(u) is identical 
to the scaled capture- number distribution C{u). Finally, 
by slightly modifying our analytical expression (|14|l for 
the scaled CND and solving self-consistently, we have ob- 
tained analytic expressions for the scaled CND and cor- 
responding ISD which are in excellent agreement with 
simulations. 

We note that our analytical results also demonstrated 
that the "self-averaging property" of the capture zones 
(i.e. the average capture number Us is exactly equal to 



the local capture number evaluated at the average cap- 
ture zone of islands of size s) is exact in the asymptotic 
limit for the case of irreversible growth of point-islands 
in one-dimension. If this "self-averaging property" holds 
generally (as suggested by the good agreement obtained 
in our numerical results for finite D/ F in both one- and 
two-dimensiong^'^) then it may be possible in general to 
obtain accurate asymptotic ISDs without having to know 
the exact distribution of areas or the redistribution of ar- 
eas due to nucleation events. This represents a significant 
simplification of the calculation of ISDs since the deter- 
mination of the full distribution of capture zones is a very 
difhcult problem for which no exact solution seems to be 
possible even in the most simple casesi^ 

Finally, we consider the extension of these results to 
other cases of interest such as growth in two-dimensions. 
The present analytical work was essentially dependent 
on the asymptotically valid relation B{u) = C{u). It is 
not clear that such an explicit relation may be derived 
for the case of two-dimensional growth, and additional 
mathematical difficulties occur due to complicated ex- 
pressions for the local capture numbers. As a result, it 
is not clear if a straightforward extension of the present 
work to this case is possible, and further work is required 
to answer this question. 
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Our model of irreversible growth of point-islands in one- 
dimension may be defined as follows. Atoms are deposited 
randomly on an initially empty line of sites with a (per 
site) deposition rate F and difi'use with a diffusion con- 
stant D. When a monomer moves onto a site occupied by 
another monomer, a dimer island is nucleated at that site. 
Similarly, a monomer moving onto a site occupied by an 
island is absorbed and the island size s increases by one. 
All islands of size s > 2 are assumed to be stable and 
immobile. 

In the asymptotic limit a vanishingly small fraction of is- 
lands will have nucleated at a coverage 9y which does not 
satisfy 6y » 9x and therefore this assumption holds for 
essentially all islands. 

For example, Eq. 1151 with n = 3 leads to a distribution 
with a slightly higher peak which is just slightly shifted to 
the right and with a tail which is just slightly lower than 
in simulations. 



